the block lsmr algorithm for solving linear systems with multiple right-hand sides

Authors

maryam mojarrab

faezeh toutounian

abstract

lsmr (least squares minimal residual) is an iterative method for the solution of the linear system of equations and leastsquares problems. this paper presents a block version of the lsmr algorithm for solving linear systems with multiple right-hand sides. the new algorithm is based on the block bidiagonalization and derived by minimizing the frobenius norm of the resid ual matrix of normal equations. in addition, the convergence of the proposed algorithm is discussed. in practice, it is also observed that the frobenius norm of the residual matrix decreases monotonically. finally, numerical experiments from real applications are employed to verify the effectiveness of the presented method.

Upgrade to premium to download articles

Sign up to access the full text

Already have an account?login

similar resources

the block lsmr method: a novel efficient algorithm for solving non-symmetric linear systems with multiple right-hand sides

it is well known that if the coefficient matrix in a linear system is large and sparse or sometimes not readily available, then iterative solvers may become the only choice. the block solvers are an attractive class of iterative solvers for solving linear systems with multiple right-hand sides. in general, the block solvers are more suitable for dense systems with preconditioner. in this paper,...

full text

The block least squares method for solving nonsymmetric linear systems with multiple right-hand sides

In this paper, we present the block least squares method for solving nonsymmetric linear systems with multiple righthand sides. This method is based on the block bidiagonalization. We first derive two algorithms by using two different convergence criteria. The first one is based on independently minimizing the 2-norm of each column of the residual matrix and the second approach is based on mini...

full text

Block Bidiagonalization Methods for Solving Nonsymmetric Linear Systems with Multiple Right-hand Sides

Many applications require the solution of large nonsymmetric linear systems with multiple right-hand sides. Instead of applying an iterative method to each of these systems individually, it is often more eecient to use a block version of the method that generates iterates for all the systems simultaneously. In this paper, we propose block versions of Galerkin/minimal residual pair of bidiagonal...

full text

A Block-QMR Algorithm for Non-Hermitian Linear Systems With Multiple Right-Hand Sides

Many applications require the solution of multiple linear systems that have the same coeecient matrix, but diier in their right-hand sides. Instead of applying an iterative method to each of these systems individually, it is more eecient to employ a block version of the method that generates iterates for all the systems simultaneously. In this paper, we propose a block version of Freund and Nac...

full text

New variants of the global Krylov type methods for linear systems with multiple right-hand sides arising in elliptic PDEs

In this paper, we present new variants of global bi-conjugate gradient (Gl-BiCG) and global bi-conjugate residual (Gl-BiCR) methods for solving nonsymmetric linear systems with multiple right-hand sides. These methods are based on global oblique projections of the initial residual onto a matrix Krylov subspace. It is shown that these new algorithms converge faster and more smoothly than the Gl-...

full text

Product Hybrid Block GMRES for Nonsymmetrical Linear Systems with Multiple Right-hand Sides

Recently, the complementary behavior of restarted GMRES has been studied. We observed that successive cycles of restarted block BGMRES (BGMRES(m,s)) can also complement one another harmoniously in reducing the iterative residual. In the present paper, this characterization of BGMRES(m,s) is exploited to form a hybrid block iterative scheme. In particular, a product hybrid block GMRES algorithm ...

full text

My Resources

Save resource for easier access later


Journal title:
iranian journal of numerical analysis and optimization

جلد ۵، شماره ۲، صفحات ۱۱-۰

Keywords

Hosted on Doprax cloud platform doprax.com

copyright © 2015-2023